Number distributions for fermions and fermionized bosons in periodic potentials 
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We compute the spatial population statistics for one-dimensional fermi-gases and for bose-gases 
with hard core repulsions in periodic potentials. We show how the statistics depend on the atomic 
density in the ground state of the system, and we present calculations for the dynamical turn-on of 
the potential. 



PACS numbers: 03.75.Fi, 42.50.Ct 
INTRODUCTION 

Since the first achievement of Bose-Einstein condensa- 
tion (EEC) of alkali gases in 1995, atomic quantum gas 
research has continued to provide important tests of fun- 
damental physics. One very ambitious current goal for 
ultra-cold atom research involves cold fermionic atoms, 
where association of the fermionic atoms into bose con- 
densed diatomic molecules has just been reported by sev- 
eral groups 1], and where current efforts aim at the ob- 
servation of BCS superfluidity of the gas |3l- Another 
goal is to push the study of cold bosons far beyond the 
mean field regime, as exemplified already when bosons 
with a repulsive interaction are exposed to a periodic po- 
tential giving rise to the transition to a Mott insulator 
[3, IJI , and in the case of a very dilute gas in one dimen- 
sion, where the repulsive interaction effectively presents 
a hard core exclusion of overlapping atomic population 
and leads to the To nks- Girardeau regime [3. 

The zero temperature state of the Tonks-Girardeau 
bose-gas (TG-gas) is described by a wave function which 
is obtained mathematically from the Slater determinant 
describing a spin-less non-interacting fermi-gas. This has 
the consequence that these two systems have many com- 
mon properties, and in the present paper we wish to 
determine one of these properties, namely the fluctu- 
ations in the number of atoms populating a finite re- 
gion of space. The population statistics in weakly in- 
teracting and in ideal dilute Bose gases was studied in 
|6(, J^ , which revealed a pair-correlation mechanism for fi- 
nite numbers of non-condensed particles. Measurements 
of fluctuations represent an interesting approach to the 
analysis and understanding of physical systems which, 
in analogy with pioneering experiments in quantum op- 
tics, extracts its information from the noise rather than 
from mean values in measurements. In other many-body 
problems, noise measurements have already been used to 
demonstrate the fractional quantum Hall effect, and the 
Hanbury-Brown and Twiss correlations in high energy 
nuclear collision experiments where they may constitute 
probes of quark-gluon plasma effects. 

We note that a TG-gas recently was produced and ob- 



served in experiments, where a periodic potential was 
used to alter the effective mass and hence the importance 
of the interaction energy relative to the kinetic energy of 
the particles |3]- The analysis in that paper was based 
on measurements of the momentum distribution. In this 
paper we develop the theory for another observable of 
this system, namely the counting noise in real space. We 
shall also subject our atoms to periodic potentials and 
investigate how this alters the counting statistics. 



SPATIAL POPULATION STATISTICS FOR 
FERMIONS AND TG-BOSONS 

A pure state of N non-interacting spin-polarized 
fermionic particles (i.e., all populating the same internal 
spin state) is described by an anti-symmetrized product 
state (Slater determinant) of single-particle spatial wave 
functions 



^^(a^i, •■■, xn) oc det{ui{xj,t)}, 



(1) 



where the N orthogonal single-particle states Ui{x,t) are 
all solutions of the single-particle Schrodinger equation of 
the atoms. The stationary ground state is thus obtained 
from the N lowest single-particle energy eigenstates. 

A pure state \E'_B(a^i, ■■■, a^Af) of iV impenetrable bosons 
must be symmetric under exchange of the atomic coor- 
dinates, and it must satisfy the condition that it van- 
ishes if any two particle coordinates coincide: ^s = if 
Xi = Xj for any i ^ j- In the limit where the hard-core 
interaction between the boson has zero range, a state 
vector "^sixi, ...,xn) obeying these constraints is for- 
mally obtained by introducing the antisymmetric func- 
tion A{xi, ...,xn) = T^i<i<j<NSga{xj — Xi), and simply 
defining, 



*B(a;i, ...,xn) = A{xi, ...,xn)'^f{xi, -..^xn)- (2) 

This 'fermionization' ansatz is valid both for stationary 
states (e.g., the ground state of the system), and for time 
dependent problems. 



When the state vector is known, the probability that 
precisely n particles are found in an interval of space /, 
is obtained by an integration over all configurations con- 
tributing precisely to this possibility, i.e., n coordinates 
should be within the interval / and the remaining N — n 
coordinates should be outside that interval: 



J Xi-^,...Xi^^I <^ Xj-^ ,...Xj, 

X |*F/B(a;i,...2;Ar)p. 



dxi . . .dx 



N 



ii 



(3) 



The sum extends over all combinations 1 < ii < 12- ■■ < 
in < N, 1 < ji < ... < JN-n < N with ykj : ik ^ ju 
and although the different contributions to this sum all 
yield the same value, one still has to carry out an A'^- 
dimensional integral of a Slater determinant with iV! 
terms to obtain the distribution. We have written the 
expression Q here because it shows explicitly that non- 
interacting fermions and hard-core interacting bosons 
have identical number distribution on a given position 
range. By writing the counting distribution in terms of 
atomic creation and annihilation operators, we shall ob- 
tain a simpler expression for the counting statistics. 

The treatment below will address the fermionic system 
and it will make use of operators with fermionic anti- 
commutator properties, but the resulting number distri- 
butions on spatial intervals are applicable to both non- 
interacting fermions and impenetrable bosons. Note that 
the state vectors in momentum space are not linked as 
easily as the expressions (I1I2II in position space: the free 
fermi gas at zero temperature has a flat momentum dis- 
tribution up to the fermi momentum, whereas the Tonks 
gas has a characteristic peaked distribution around zero 
momentum. In (^ no single-particle state is populated 
by more than a single atom, whereas the Tonks gas has 
a single state populated by of the order vTV atoms. 

Let a„ denote annihilation operators for fermions in 
a state with single-particle wave function u„(a;). The 
connection between these operators and the position de- 
pendent particle annihilation operator 7/1(2;) is 



Vi(x) 






Un{x)a„ 



dxu*Jx)il}{x) (4) 



where we make use of both orthogonality and complete- 
ness of the basis of eigenstates Un{x). 

In the calculations below, we shall assume that the 
single-particle states n = 1, . . . ,N are all occupied. The 
probability to detect a single atom at the location xi 
is the expectation value of the local density operator 
ip''{xi)'4'{xi), and we can evaluate the average popula- 
tion of a position interval / as the mean value of the 



number operator fij — J , dxiilj''{xi)ip{xi), 



{ni) = 



N 

E 



dxi\Un'{xi)\" 



(5) 
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Making use of the anti-commutator ?/;(a:)-0'l'(a;') _(_ 
ijj''{x')ip{x) = S{x — x') we can also address higher or- 
der moments of this population. 



N 



E / dxi\Un'{xi)\'^ 

N n'-l „ „ 

E E ' ^^1 ' ^^2 



n' — 1 n" — 1 



xiEl 



X2el 



x\Un'ixi)Un"{x2) - Un' {x2)Un" {xi)\'^ . (6) 

If we know the single-particle states occupied by the 
system, we can compute the mean population and the 
variance in the interval / explicitly. For example, in 
the lowest energy state of free atoms, where all the 
eigenstates are plane waves, we obtain for the number 
of fcrmions/TG-bosons on a ID interval of length L: 
{•ill) — pL, where p is the mean density of the gas, and 
we straightforwardly obtain a variance that scales with 
the logarithm of the mean occupancy of the interval. 
For an infinitely extended system, one can also obtain 
this result from the Fermi-Dirac momentum distribu- 
tion function, which is the Fourier transform of the spa- 
tial density-density correlation function j9(], Var(n7) — 
(1 -I- ln(27rpi) + 7)/7r2, with Euler's 7 = 0.57721... . In 
quantum optics language, the fermions and TG-bosons 
are strongly anti-bunched. 



Full number distribution 

By use of a technique developed by Levitov [iflj (for 
a simpler derivation of this single aspect of Levitov's re- 
sults, see the presentations by Klich llj and by Castin 
1^ Il3l)i it is possible to obtain the full number dis- 
tribution p{n) for the number of atoms in the inter- 
val /, from which all moments of the distribution are 
readily calculated. First, we observe that the operator 
■^ L dd ex'p{i9{ni — n)) is the projection operator onto 
a subspace where the operator fij has the eigenvalue n, 
see also |3- It is therefore useful to introduce the char- 
acteristic function 



Fi9) = (exp{i9ni)), 



(7) 



from which the number distribution is obtained as a sim- 
ple Fourier transform. The expectation value in Q is 
taken in the state with N occupied single-particle wave 



functions, and the trick is to realize that this state is de- 
scribed by a density matrix that can formally be written 
as a thermal state a = exp{—(3H)/Z where f3 = l/ksT, 
and where H — J2 ^nO-li'^n is the many-body (second- 
quantized) Hamiltonian corresponding to a one-body 
Hamiltonian H (for clarity we write operators acting on 
the one-particle Hilbert space without hats), which is di- 
agonal in the basis of states u„ (x) . The expectation value 
in (0 is thus the trace of a product of two exponen- 
tials of quadratic forms in field creation and annihilation 
operators. In general, if operators ^4, B and C on the 
one-particle Hilbert space obey exp(A) exp(B) = exp(C), 
then also the relation exp(A) exp(_B) = exp(C') is valid 
for the second-quantized expressions X = ^ Xnm0.nO,m, 
X ^ A, B, C. Apart from the factor 1/Z, our character- 
istic function is thus the trace of an operator exp(C') — 
exp(^^ Xiblbi), where A^ are the eigenvalues of the ma- 
trix C, and bi are the corresponding linear combinations 
of the operators a„. Since the &|&i operators are Fermi 
number operators, their eigenvalues are restricted to the 
values zero and unity, and the exponential exp(Ai6j6i) 
yields the values 1 and exp(Ai). Hence we can write 



Tr[e^e^] = Tr[e'^] 

= Tr[exp(^A,6|6,)] 



(8) 



= Tr[H,cxp(A,6|6,)] 



H.(l- 

det[l 

det[l 



e^e^l 



expressing the trace of exponentials of second-quantized 
operators in terms of a determinant involving exponen- 
tials of one particle operators [11|. We now insert the 
relevant matrices, 



F{e) = Tr[e""'^^^e-^^] = ^ det[l + e^'^^'e'^^ (9) 

Z/ Zj 

where Pi — ^ -,\x){x\ is a projection operator 
(p2 ^ p^)^ and thus obeys e'**^^ = 1 -h (e*^ - 
1)P/. Introducing the population of the single-particle 
states, 7r„, the characteristic function then rewrites 
F{Q) = det[l + {e'" - 1)P/ Y^^ 7r„|n)(n|], and if we repre- 
sent this operator as a matrix in a discrete position rep- 
resentation with grid spacing Ax, so that the projection 
Pj effectively extracts only a sub- matrix with x,x' € /, 
we get 



Schrodinger equation, this expression allows us to de- 
termine the population statistics. As stated above, the 
expression is valid both for non-interacting fermions and 
for TG-bosons at zero temperature, where most of our 
calculations are performed. In this case, 7r„ = 1 for all 
occupied states. We shall also study how the statistics 
depends on temperature. However, because the Bose- 
Fermi mapping only holds for T = 0, these results only 
apply to fermions. 

Since the determinant is formally a linear combina- 
tion of products of Nx matrix elements, where N^ is the 
number of x-values used to represent the interval /, F{6) 
becomes a polynomial in e*^ of order N^. This number 
should be large enough to resolve the spatial structure of 
the wave functions populated, and we see that the Fourier 
transform assigns vanishing probability to any number of 
particles larger than N^, which is meaningful when we 
represent the interval by N^ localized basis states. For 
weakly interacting bosons, their commutators lead to a 
generating function which has the e'^ dependence in the 
denominator |3, lll| . and hence any high occupancy is 
in principle possible. In our numerical calculations, we 
shall of course check the convergence of our results when 
larger and larger iV^; values are applied. 

For comparison, we note that the number distribution 
of a pure condensate of non-interacting bosons is bino- 
mial, i.e. 



p{n) = 



Nl 



n!(iV 



-y\i-py 



(11) 



where p — J , da;|uoP is the probability for a particular 
particle to be in /, as given by the macroscopically popu- 
lated single-particle wave function uq in the condensate. 
This result holds irrespective of the external potential, 
and in particular, we note that the number distribution 
for a single lattice site is independent of the amplitude 
of the periodic potential. This ideal bosonic behavior 
does not apply when the bosons interact. In the case 
of repulsive interactions the condensate exhibits a phase 
transition to a Mott-insulator with a specific number of 
particles pr. lattice site when the lattice amplitude is 
raised above a certain value yi Ijl . 



PARTICLES IN A PERIODIC POTENTIAL: 
NUMERICAL RESULTS 

We shall now use the theory for the spatial statistics 
to calculate the number distribution of fermions and TG- 
bosons in a single well of a periodic potential. 



Fie) = det[l -I- (e^'^ - l)AxY,T^nKix)unix')]. (10) 



Statistics of the ground state 



Knowing the single-particle states, either analytically, by 
diagonalization or by propagation of the time-dependent 



We first investigate the ground state properties of the 
system, i.e. we assume that the N particles occupy the N 



a) 







FIG. 1; Atom number distributions for fermions and TG 
bosons in a single potential well as a function of the lattice 
amplitude Vo. a) mean occupation n — 5 and b) mean occu- 
pation n = 4.4. The distributions were calculated assuming 
r = and that the periodic potential was turned on adiabat- 
ically. The amplitude is given in units of the lattice energy 
Eiattice defined in the text. For comparison, the distributions 
for non-interacting bosons (dark grey bars) are shown at the 
same mean occupancies. 



single-particle states with lowest energy. To find the wave 
functions Un{x) of these states required to obtain the 
number statistics by Fourier transformation of l|10() . we 
perform a band structure calculation to solve the single- 
particle Schrodinger equation 

--Z -ri- + -^ n)(l + COS[Kx\)Un^q = hhJn,qUn^q ■ (12) 

Here, m is the atomic mass, and k = 47r/A, where A is 
the wavelength of the optical lattice beams used to cre- 
ate the periodic potential. Vq is the amplitude of the 
periodic potential, which is proportional to the intensity 
of the lattice beams and therefore adjustable experimen- 
tally. The natural unit for the amplitude is the lattice 
energy Eiatuce — 2/i^/(mA^), which differs from the re- 
coil energy associated with the optical lattice beams by 
a factor of four. Due to the periodicity of the Hamilto- 



nian, the wave functions take the Bloch form Un,q{x) = 
exp{iqx)4>n,qix), where the integer n is the band index 
and q is the crystal momentum, which under the assump- 
tion of periodic boundary conditions over n^ periods, 
can take riw discrete values usually taken to be in the 



first Brillouin zone: q 
El. The function d 



„/2 



e{i,2,. 



} 



',^q{x) is periodic with a spatial 
period A/2. We expand it on a plane-wave basis set 
(bn,,{x) = E'Sr^'c„,„,((7)exp[*f x(n, - n™- - 1)], 



and solve the Schrodinger equation by finding the coef- 
ficients Cn,nb (9) of the eigenvectors of the Hamiltonian 
in this basis. Once the eigenvectors are known for suffi- 
ciently many q, and as a function of Vq, the effect of the 
periodic potential on the number statistics can be calcu- 
lated by insertion of the wave functions into 1)10(1 and by 
calculating the Fourier transform of F{9). 



Figure ^ shows the number distribution for fermions 
and TG-bosons at a lattice site as a function of Vq for 
two selected mean occupancies n — 5 and n = 4.4. 
Also shown in the figure are the distributions of non- 
interacting bosons at the same mean occupancy (dark 
grey bars). For free particles (Vq = 0), the number dis- 
tribution of fermions and TG-bosons is already much nar- 
rower than that of non-interacting bosons. Moreover, the 
distribution of TG-bosons and non-interacting fermions 
changes with Vq, in contrast to non- interacting bosons. 
In the case of commensurate filling (n integer), the dis- 
tribution becomes P{n) = SnM for large Vq, and conse- 
quently the variance of the atom number goes to zero. 
This is exactly what one would expect for fermions be- 
cause the Pauli principle makes it energetically favor- 
able to distribute the particles evenly over the lattice 
sites. When the filling is incommensurate (n = riint + e, 
where riint is the integer-part of n) , the distribution be- 
comes bimodal at large Vq, with P{nint) = (1 — e) and 
P{nint + 1) = e. Hence, the variance converges towards 
the value Var/ — e(l — e), which is symmetric around 
e = 0.5 and independent of riint- The dependence of 
Var/ on Vq is shown in Fig.|21for a selection of mean oc- 
cupancies in the range 3 < n < 4. In addition to the sym- 
metry of the variance around e = 0.5 mentioned above, 
the figure shows that the lattice amplitude required for 
the probability distribution to converge increases with 
the mean occupation. To address this issue quantita- 
tively, we calculated the amplitude V™^ required for 
the variance to be reduced to 1% of the variance for free 
particles as a function of n for integer n. Our findings 
are listed in Tabled which also compares V^^''^ to the 
mean energy of the highest occupied band {E)band at a 
lattice amplitude Vq™"^"^^". As one might expect, the lat- 
tice amplitude has to be slightly larger than the mean 
energy of the highest occupied band in order to reduce 
the variance to the few-% level. 



0.5 



0.4 



0.3 



0.2 



0.1 



^v 






\\\X^^V 






\\\3.3V^^^ 
\ \ 3.2A \ 


^^^i6r===^ 


3.5 


:i..,I^^3.7- 




^^ 


- 


\ •^•'' \ 


\--_ 3.8~_____ 




""•v 




3.0 


\___ 3.9._ 




\ 




v^ 


4.0 





at that particular time: 



12 

^0 L^isWcs J 



16 



20 



FIG. 2: Variance of the number distribution for fermions and 
TG-bosons at a single lattice site as a function of lattice am- 
plitude Vb for various mean occupancies n = riint + e. The 
variance converges towards e(l — e) at large amplitudes. The 
amplitude required for convergence increases with n. The 
calculations assume that the periodic potential is raised adi- 
abatically and that T = 0. 



'. 12 3 4 5 6 7 8 

[Ei^ttice\ 2.5 5.9 10.7 16.8 24.1 32.8 42.7 53.8 
{S>band [Eiatuce] 0.7 3.3 7.2 12.5 19.1 26.9 36.1 46.5 



T/var 



= 1% 



TABLE I: The row VJ""'"^''" Usts the lattice amplitude re- 
quired to reduce the variance of the atom number distribution 
at a lattice site to 1% of the variance for free particles as a 
function of the mean occupancy n. The row (i5)band contains 
the mean energy of the highest occupied band at a lattice 
amplitude of Vq'^'' . The energy unit is Eiatuce- 



Effect of non-adiabatic turn-on 

So far we have only considered the ground state of 
the system, i.e. we have assumed that the periodic po- 
tential is turned on adiabatically, so that only the N 
single-particle states with lowest energy are occupied at 
all times. In reality, the potential will be turned on using 
a specific ramp Vo{t), and if the ramp is too fast, atoms 
will be driven into initially unoccupied bands, which will 
lead to an increase in the atom number fluctuations at 
a lattice site. We shall now analyze how the lattice 
turn-on affects the state and the spatial statistics of the 
many-body system. To this end, we calculate the time- 
evolution of the state of the system for various ramps by 
solving the time-dependent Schrodinger equation for the 
single-particle states. 

Our procedure is as follows: At each time t, we expand 
the single-particle states (/>$^f^(i) on the Bloch eigenfunc- 

tions </>„'„/ [Vo(i)] at the amplitude Vo[t] set by the ramp 



€':{t) 



/ . "n,q;n',q'(Pn',q' C " " 



(13) 



where hujn'^q'(t') is the energy of the state (j^ni „,[Vo{t')]. 

This expansion is inserted into the time-dependent 
Schrodinger equation which results in a set of first-order 
differential equations for the time dependence of the co- 
efficients bn^q-n'y. 

To solve these differential equations, we need to know 
the matrix elements (</>„/ „/ 1^^„// „//) and the band struc- 
ture LUn'.q' as a function of Vo(0- We therefore start our 
calculations by solving the static Schrodinger equation 
(|12|l for a sufficiently large number of Vq 's ranging from 
zero to the final Vq of the ramp. We then calculate the 
matrix elements: 



n ,q 



d 



,(0) 






dt 



"6 = 1 



M' 



dVo 



(14) 

Here, c„'^„j,(g') are the plane wave expansion coefficients 
introduced previously. Note that the matrix elements 
are diagonal in the crystal momentum quantum num- 
ber and, consequently, the dynamics can be solved for 
each q independently. Once the matrix-elements and en- 
ergies of the differential equations have been calculated, 
the coefficients hn,q-^n',q' are propagated in time. We as- 
sume that the system is in the ground state initially, i.e. 
hn,q;n',q'{t^^) = 5n,n'5q,q' for {ti, g} Corresponding to 
the N states with lowest energy and bn,q-n\q'{t = 0) = 
for all other states. 

To quantify the degree of excitation caused by the 
turn-on, the time-evolved single-particle states are used 
to calculate the population of those Bloch states, which 
are not populated initially. To make this number in- 
dependent of the number of periods riw for the peri- 
odic boundary condition, we shall discuss the degree 
of excitation in terms of the number of excited atoms 
pr. period of the periodic potential A^oxc- To investi- 
gate how the excitations affect the atom number statis- 
tics, we also insert the time-evolved single-particle states 
into Hl()() . There are infinitely many ways to turn on 
the optical lattice. We shall limit our discussion to two 
cases: linear and exponential ramps. The linear ramps 
start at Vq = and increase linearly to the final value 
Vo. final = QOEiatuce over a time interval of length iramp- 
In the experimentally relevant case of '^''K fermions and 
an optical lattice wavelength of A = 850 nm, we shall 
study ramp times up to 80 h/Eiatuce ~ 461 /iS. Our 
exponential ramps have the time-dependence Vo(t) = 

Vb,/mai[exp(5t/tramp)-l]/[exp(5)-l], i G [0, tramp]- This 

ramp is very similar to the ones used in the recent exper- 
imental realization of a TG-gas |3] ■ 

To study the cross-over from nonadiabatic to adiabatic 
turn-on and the effect on the atom number statistics for a 
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FIG. 3: Number of excited atoms pr. lattice site as a function 
of tiie ramp time and of the mean number of atoms pr lattice 
site (for commensurate filling). The solid (dashed) curves are 
the results for linear (exponential) time dependencies. For 
clarity, the curves for different n are offset vertically with a 
spacing of AA'oxc = 0.1. Thus, the spacing between horizon- 
tal lines defines the scale for the vertical axis. The duration 
of the ramp is given in units of h/Eiattice, where Eiatuce is 
the lattice energy defined previously. For ''"K and an optical 
lattice wavelength of A = 850 nm, Ti/Eiatuce ~ 5.8/iS 



single lattice site, we calculated the time evolution for lin- 
ear and exponential ramps with durations 5h/ Eiatuce < 
iramp < ^GH/ Eiattice and for vailous mean occupancies n. 
The final amplitude was constant Vojmai = (^QEiatuce for 
all calculations presented here. Fig. |31 shows the num- 
ber of excited atoms pr lattice site (Nexc) as a func- 
tion of tramp for integer mean occupancies in the range 
1 < n < 8. The solid curves are the results for lin- 
ear ramps and the dashed curves represent exponential 
ramps. The curves corresponding to different n are offset 
vertically from each other by ATVexc = 0.1, and for each 
n, the thin horizontal solid line corresponds to A^cxc = 0. 
Thus, the spacing between the horizontal lines defines the 
vertical scale of the figure. Starting in the lower right cor- 
ner of the figure, we see that for large mean occupancies 
fi G 7,8, and ramp times larger than ~ AQh/ Eiatuce, the 
number of excited atoms pr lattice site is at or below the 
few-% level. Thus, the turn-on is essentially adiabatic. 
Moreover, linear ramps lead to less excitation than the 
corresponding exponential ramps. This also holds for 
faster ramps, where the degree of excitation increases. 
When the mean number of atoms pr site is decreased to 



4 < n < 6, the excitations caused by the two types of 
ramps are quite similar. For small mean occupancies, 
n < 3, exponential ramps are superior because they rise 
more slowly at early times where the ramp surpasses the 
energy of the occupied bands. 

Our adiabatic calculations show that (cf. Fig. |2Jl par- 
tially occupied bands (incommensurate filling) lead to 
fluctuations in the number of atoms pr lattice site even 
at large Vq- Non- adiabatic turn-on of the periodic po- 
tential transfers a fraction of the atoms from the highest 
occupied to the lowest unoccupied band, and therefore re- 
sults in (at least) two partially filled bands, which leads to 
fluctuations in the site occupancy also in the case of com- 
mensurate filling. Using the time-evolved single-particle 
states and our theory for the spatial statistics, we find 
that the probability distribution at large Vq for integer n 
broadens compared to the adiabatic monomodal distri- 
bution p{n) = Sn.n- Except for (a) very fast exponential 
ramps (tramp < ^Q^/ Eiattice) and (b) ramp durations of 
tramp < ^^fi/ Eiattice and Small mean occupancies n < 2 
for linear ramps, the distribution is to a good approxima- 
tion given by {p{fi— 1) — rj,p{n) = 1 — 2ri,p{n+ 1) = 77}, 
where 77 = iVexc- This distribution has a variance of 2?^, 
and the variance can therefore be obtained from Fig. |31 
simply by scaling iVoxc by a factor of two. 



Effect of non-zero temperature 

All the results presented so far correspond to zero tem- 
perature. In this last section we shall discuss how non- 
zero temperature affects the spatial statistics of fermions 
in periodic potentials. The bose-fermi mapping |(2l, 
which is the reason that the spatial statistics of TG- 
bosons is identical to that of fermions, does not apply 
at non-zero temperature. Our theory and the results 
presented in this section does therefore not apply to TG- 
bosons. 

It is straightforward to generalize our numerical simu- 
lations to non-zero temperature. The only difference lies 
in the occupation probabilities of the states 7r„ in H10|l . 
which for non-zero T are given by the Fermi-Dirac distri- 
bution: TTn^q = (exp[/?(?ki;„^g — fi)] + 1)~^, where /z is the 
chemical potential defined by the constraint that the sum 
of TTn.q ovcr all states equals the appropriate total num- 
ber of particles. We assume that the system is in thermal 
equilibrium with a temperature T and at the amplitude 
Vb, i.e. our results correspond to the case were the sys- 
tem has equilibrated at Vq. Alternatively one might have 
assumed that the system equilibrates so slowly that it re- 
tains the distribution iTn.q corresponding to free particles 
{Vq = 0), and then calculate the statistics using these 
occupation probabilities. We shall only consider the first 
case here. 

Assuming that the periodic potential has been raised 
adiabatically, there are two contributions to the fluctu- 



ations in the number of particles at a lattice site. In 
addition to the fluctuations due to the fermionic nature 
of the particles, which is present at T = and was dis- 
cussed above, thermal excitation lead also lead to fluc- 
tuations. This can be seen in Fig. 01 which shows the 
variance as a function of Vq for two mean occupancies 
{n = {2, 4}) and as a function of temperature. The tem- 
perature is given in units of Eiatuce/kB- In the case of 
'^"K and A = 850 nm, the four temperatures are nK, 
120 nK, 240 nK, and 360 nK. According to the figure, 
the variance increases with temperature for a given Vq 
due to thermal excitations. For small Vq, the increase in 
variance due to thermal excitations decreases when par- 
ticles are added to the system. The reason is that the 
mean spacing between bands increases with increasing 
band index. Since the Fermi-Dirac distribution changes 
from unity to zero over an energy range kBT, the increase 
in energy level spacing leads to less thermal excitation. 
Fig. 01 also shows that the variance, for the case of com- 
mensurate filling, goes to zero when the temperature is 
non-zero. The amplitude Vq required to obtain a certain 
variance increases with T. The reason for this behav- 
ior is also easily understood. At large Vq, the particles 
are localized near the periodic potential minima, and the 
potential can be expanded to second order around the 
minima. In this regime the bands a equally spaced with 
an energy spacing proportional to VMi- Thus, the en- 
ergy spacing between bands increases with Vq and, con- 
sequently, the thermal fluctuations vanish for large Vq. 

The thermal contribution to the atom number distri- 
bution can be studied separately by choosing Vq such 
that the variance at T = is much less than the vari- 
ance at the temperature of interest. In this regime, 
the atom number distribution for commensurate filling 
is {p(n — 1) = r],p{n) = 1 — 2r],p{n + 1) = ij}, like in the 
case of excitations caused by non-adiabatic turn-on. In 
the thermal case, i] is equal to J2n n ''^n.q/'riw, where the 
sum extends over all the states not occupied at T = 0. 



CONCLUSION 

We have presented a theory for the spatial statistics 
of fermions and Tonks-Girardeau bosons in one dimen- 
sion, and applied the theory to the case of fermions and 
TG-bosons in periodic potentials. The number distribu- 
tion of fermions and TG-bosons occupying a single lattice 
site is significantly narrower than that of non-interacting 
bosons. Moreover, the distribution depends on the am- 
plitude of the periodic potential and becomes bimodal 
(monomodal for commensurate filling) at amplitudes on 
the order of the energy of the highest occupied band. 
This behavior is completely different from that of non- 
interacting bosons, which is independent of the potential 
depth. 

We also studied the dynamics associated with the turn- 
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FIG. 4: Variance of the femiion number distribution as 
a function of Vo for n = {2, 4} and temperatures T = 

{0,0m,0.W,0.24:}ElatUce/kB. 



on of the periodic potential and the effect of excita- 
tions on the number distribution of a lattice site. Non- 
adiabatic turn-on of the potential leads to partially filled 
bands and fluctuations in the site occupancy even in the 
case of commensurate filling. For commensurate filling 
and most ramps of practical interest, the probability of 
having one atom more or less than the mean is equal to 
the average number of excited atoms pr lattice site. The 
variance is well approximated by twice this number. 

Finally, we studied how thermal excitations affect the 
statistics of fermions in periodic potentials. We found 
that the fluctuations in the number of atoms at a lattice 
site increase with temperature. The variance of the num- 
ber distribution converges towards the variance at T = 
when the potential depth is increased, but the potential 
depth required to obtain a certain variance increases with 
temperature. For commensurate filling and lattice depths 
large enough that the fluctuations at T = are vanish- 
ingly small, thermal excitations lead to a symmetric atom 
number distribution. Like in the case of non-adiabatic 
turn-on, the probability of having one atom more or less 
than the mean is to a good approximation equal to the 
mean number of excited atoms pr lattice site. 

Yvan Castin is acknowledged for stimulating discus- 
sions and for referring us to the work on the general pop- 
ulation statistics. 
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